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Abstract. Observations of distances to Type-la supernovae can be explained 
by cosmological models that include either a gigaparsec-scale void, or a cosmic 
(— I ' flow, without the need for Dark Energy. Instead of invoking dark energy, these 

Q_c inhomogeneous models instead violate the Copernican Principle, we show that current 

q \ cosmological observations (Supernovae, Baryon Acoustic Oscillations and estimates 

H ' °f the Hubble parameters based on the age of the oldest stars) are not able to 

rule out inhomogeneous anti-Copernican models. The next generation of surveys 
for baryonic acoustic oscillations will be sufficiently precise to either validate the 
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Copernican Principle or determine the existence of a local Gpc scale inhomogeneity. 
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tyQ ' 1. Introduction 

O" 

The Copernican Principle states that we do not occupy any special place in the Universe. 
The Friedmann-Lemaitre-Robertson- Walker (FLRW) models are built on this principle, 
and provide a remarkably precise description of cosmological observations [H, El E]. 
One could therefore think of this as an indirect demonstration of the Copernican 
Principle. However, homogeneous and isotropic FLRW models are not the only solutions 
of Einsteins equations which are able to fit cosmological observations. In particular, 
a number of inhomogeneous models have been proposed, each of which are able to 
describe the evolution of distance with redshift as measured via type-la supernova (Snla) 
without need of a cosmological constant [U [5j EJ U\ E] (see [9] for a review) . Moreover, 
these models may be constructed in such a way that they describe the details of the 
cosmic microwave background power spectrum (CMB) [5], [6]. By eshuing the constraint 
of inhomogeneity, these inhomogeneous models violate the Copernican Principle and 
suggest that we live near the center of a large (of several gigaparsecs diameter) highly 
isotropic void. 

Inhomogeneous models are able to fit a variety of sets of cosmological observations 
without containing a cosmological constant because the lack of homogeneity offers 
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a great degree of flexibility [TO]. For example, since the last scattering surface is 
separated from regions where supernova are observed by great distances, the property 
of inhomogeneity allows a model to be constructed which provides different physical 
densities in the regions from which these two sets of observational data are drawn. Thus, 
what is required to constrain inhomogeneous solutions to the apparent acceleration of the 
Universe includes several sets of data that measure a range of observables at comparable 
redshifts. In this paper we confront two general classes of inhomogeneous, dark energy 
free models with observations of Snla, with measurements of the baryon acoustic 
oscillations (BAO), and with the variation of the Hubble parameter with redshift. All 
these observational data are drawn from z < 1.8. We begin by describing the Lamaitre- 
Tolman model (the simplest spherically symmetric inhomogeneous generalization of 
FLRW model), and then construct a range of 2-parameter models to describe a 
local large scale inhomogeneity. We then discuss constraints on the scale of these 
inhomogeneities based on existing and forthcoming data. 

2. The Lemaitre-Tolman model 

The Lemaitre-Tolman (LT) model [TT] is a spherically symmetric, pressure free, 
irrotational solution of the Einstein equations. Its metric has the form 

ds 2 = c 2 dt 2 - ^^ dr 2 - R\t, r)dfi 2 , (1) 

where dQ 2 = d9 2 + sin^d^ 2 . Because of the signature (+, — , — , — ), the E(r) function 
must obey E{r) > —1/2. Here a prime (') denotes d r . 
The Einstein equations reduce to the following two 

, , 2 2M'(r) 

Kp ^ t)0 = R 2 (r,t)R'(r,ty (2) 

^(r.^^CD + g + lA^Cr.t), (3) 

where M{r) is another arbitrary function and k = %nG/c A . Here a dot Q denotes d t . 
When R' — and M' ^ 0, the density becomes infinite. This happens at shell crossings, 
and is an additional singularity to the Big Bang that occurs at R = 0, M' ^ 0. By 
setting the initial conditions appropriately the shell crossing singularity can be avoided 
(see [12] for detail discussion). 

Equation J3]) can be solved by simple integration: 

ft 

c[t-t B (r)], (4) 



2E+ 2M + 1 A ^2 



where ts appears as an integration constant and is an arbitrary function of r. This 
means that the big bang is not a single event as in the FLRW models, but occurs at 
different times at different distances from the origin. To define a particular LT model 
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two functions must be specified as initial conditions. For completeness we note that 
in the FLRW limit R(r,t) = ra{t) [where a(t) is the scale factor], M{r) = M r 3 , and 
E(r) = —k r 2 . 

3. Parametrization of Inhomogeneous Models 

In this paper we consider two families of models, in which the cosmological constant is 
set to zero (A = 0) and the observer is situated at the origin. The most popular way 
to explain Snla without A is to postulate a large scale void centered near our position 
[U El El E> IE] • The density distribution at the current instant in this first family may be 
parametrized by (see Fig. [T]) 

„2 



p{t , r) = p 



r 
1 + 5 p -5 p exp ( 



(5) 



where p — 0.3 x (3Hq) / (8irG) . In these models the big bang is assumed to occur 
simultaneously at every point (is = 0). The functions M and E are then calculated 
using eqs. ([2]) and (J3J) respectively. 

A second possible explanation for the apparent acceleration observed using Snla 
is inspired by the so called Hubble Bubble. In the Hubble Bubble phenomenon the 
expansion rate at distances beyond 100 Mpc is postulated to be slower than it is locally 
[T3] . Such a phenomenon could be described by a second family of models parametrized 
using (see Fig. [T]) 



H T {to, r ) = x = H ° 



,j*Z 



(6) 



1 - S H + S H exp I - , 

In these models density is assumed to be homogeneous at the current epoch. The 
function M can be calculated using the above relation and eq. (J3j) . As can be easily 
checked, these two families of models are regular at the origin (for origin conditions in 
the LT models see [14]). 

In each case the radial geodesies and redshift are calculated from [16] 



ln(l + z)= / 




dr— 

V1 + 2E 


--{1 + zfR. 





dt R , . „ . f , _ R 

— = . ln(l + z) = dr — (7) 

dr VI + 2E ' ' ■'■" 

and the luminosity distance is D L 

4. Cosmological Observations 

In this paper we confront the inhomogeneous models described above with three sets of 
observations. Firstly, we consider observations of Type-la supernova, which are taken 
from the Union data set [17]. In addition we also present the constraints coming from 
the Riess gold data set [TH], which until recently has been the most popular sample for 
testing inhomogeneous models. The second set of cosmological observations comprise 
measurement of the dilation scale of the BAO in the redshift space power-spectrum of 
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Figure 1. Left panel: The parameterization of the density distribution within the 
first family of model at the current instant. Right panel: The parametrization of the 
expansion rate (Ht) at the current instant within the second family of models, as well 
as the corresponding Hubble parameter Hr. The plotted profiles were obtained using 
parameters presented in Table [U 



46,748 luminous red galaxies (LRG) from the Sloan Digital Sky Survey (SDSS). The 
dilation scale is defined as 



D 



v 



n 2 



CZ 

W) 



1/3 



(8) 



H 



R 



(9) 



where Da is the comoving angular diameter distance and H(z) is the Hubble parameter 
in function of redshift. The measured value of the dilation scale at z = 0.35 is 1370 ± 
64 Mpc [19JJ. The Hubble parameter in equation (jSJ) is related to expansion along the 
radial direction, and in the LT model is given by 

B!_ 

R r 

In addition to the geometric measurements described above, the Hubble parameter 
has been estimated as a function of time based on the age of the oldest stars observed 
in galaxies at different redshifts [20]. Using eq. ((2J) it can be shown that within the 
LT model the Hubble parameter measured this way is also equal to Hr. Given a 
particular parametrization the null geodesic equations can be solved to calculate the 
luminosity distance, dilation scale and the Hubble parameter (Hr). These quantities 
can then be compared with observations to test the viability of the model by constraining 
parametrizations via a least square fit. 

| It should be noted that the value of 1370 ± 64 Mpc was obtained within the framework of the linear 
perturbations imposed on a homogeneous FLRW background. It is still an open question whether such 
analysis is appropriate or should be carried out instead in the LT background. However, we proceed 
with the value for the Dy provided by Eisenstein et al. We will come back to this issue in Sec. 16.11 
where we will try to partially estimate the errors which arise from application of different background 
models. 
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Figure 2. Likelihood profiles for Hq, a, and 5 P based on the BAO+H(z)+Union data 
sets. Left panel: The Hq likelihood profile obtained after marginalization over a and 
S p . Center panel: The a likelihood profile obtained after marginalization over Hq and 
S p . Right panel: The S p likelihood profile obtained after marginalization over a and 
Hq. The means are presented in Table [U 




Figure 3. Observational constraints on the parameters of the first family of models 
(cosmic void models), shown as contour plots of % 2 f° r the model given the data, 
after marginalizing over Hq. The confidence levels are set assuming a % 2 distribution. 
Upper Left: constraints from the BAO. Lower Left: constraints from H(z). Upper 
Center: constraints from the Union data set. Lower Center: constraints from the 
Riess gold data set. Upper Right: joint constraints from the BAO+H(z)+Union data 
sets. Lower Right: joint constraints from the BAO+H(z)+Riess data sets. 



5. Testing the Copernican Principle 

5.1. Cosmic void models 

We begin by constraining the first family of inhomogeneous models, characterized by 
t B = 0. As can be seen from equation (jSJ), the density increases from p at the origin 
to p = (1 + S p )pb at infinity. The left panel of Fig. [Q shows an example profile for 
which S p = 4.05 and a = 2.96. We note that the profile need not be extrapolated 
to infinity. Indeed, owing to the flexibility of inhomogeneous models the profile could 
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be modified arbitrarily at larger distances in order to fit other types of observations, 
including those of the CMB (see Appendix for more details). We allow the following 
three parameters to vary within specified ranges H G [64, 76] km s _1 Mpc -1 , a G [0, 6] 
Gpc, and 5 P G [0,6]. The likelihood distributions for these parameters are presented in 
Fig. El The expected (mean) values derived from the likelihood profiles are presented 
in Table [TJ The constraints on the parameters 5 P and a are shown as contour plots of 
X 2 after marginalizing over H Q in Fig. [3l As can be seen the Union data set places 
tighter constraints on the models than the Riess data set. However, it should be 
emphasized that the uncertainties in the Union data set used here do not include 
systematic errors, which may be present owning to the joining of several supernova 
surveys. Some of the systematics are associated with particular surveys, while others 
are common to all surveys (intristic variation in supernova explosion mechanism or 
their possible evolution). In the standard approach, systematic errors are estimated 
by adding a systematic component, <J sys , in quadrature to the statistical error. The 
amplitude of a sys is then evaluated so that the x 2 P er degree of freedom for the best 
fitting cosmological model is unity. However, since in this paper we aim to constrain 
inhomogeneous cosmological models, and to determine whether or not they are able to 
provide a successful fit to the cosmological observations, we have decided not to account 
for possible systematical errors. Indeed, using the procedure described above would 
allow us to fit the supernovae data with almost any model. On the other hand, we find 
that the best-fit model of the family of models studied in this section (see Table [2] for 
the exact values of model parameters) fits supernova data with x 2 — 325.89 (for 304 
degrees of freedom), which means that the model can be ruled out at only the 81.4% 
level. The overall fit to all three cosmological data sets (BAO+H(z)+Union) is 336.55 
(for 314 degrees of freedom), which means that the model can be ruled out at only the 
81.7% level. These low values of values of x 2 P er degree of freedom imply that these 
observational data cannot rule out anti-Copernican models describing a large cosmic 
depression. Within the parametrization assumed, the best fitted void has a density 
contrast of 5 P ~ 4.5 and a radius of a ~ 3.2 Gpc. 

We find that the best constraints come from the supernova data set. This is mostly 
because 304 measurements are available, compared with only 9 for H(z) and just one 
for BAO. However, as may be seen from the lower left panel of Fig. 3., the H(z) 
data prefers large voids, and so these data also contribute to constraints on parameter 
estimates. The weakest constraints come from the BAO measurement — see upper left 
panel of Fig. 3. 

5.2. Cosmic flow models 

As above we allow three parameters to vary in the range H G [64, 76] km s _1 Mpc -1 , 
a G [0,6] Gpc, and 5h G [0,0.4]. The likelihood distributions for these parameters, 
based on the BAO+H(z)+Union data sets are presented in Fig. BJ The constraints on 
the parameters 5h and a are shown as contour plots of x 2 after marginalizing over Hq 
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Figure 4. Likelihood profiles for Hq, a, and 5 P based on the BAO+H(z)+Union data 
sets. Left panel: The Hq likelihood profile obtained after marginalization over a and 
S p . Center panel: The a likelihood profile obtained after marginalization over Hq and 
S p . Right panel: The S p likelihood profile obtained after marginalization over a and 
Hq. The means are presented in Table [U 
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Figure 5. Observational constraints on the parameters of the second family of models 
(cosmic flow models), shown as contour plots of x 2 f° r the model given the data, 
after marginalizing over Hq. The confidence levels are set assuming a % 2 distribution. 
Upper Left: constraints from the BAO. Lower Left: constraints from H(z). Upper 
Center: constraints from the Union data set. Lower Center: constraints from the 
Riess gold data set. Upper Right: joint constraints from the BAO+H(z)+Union data 
sets. Lower Right: joint constraints from the BAO+H(z)+Riess data sets. 



in Fig. [5l This family of models provides an even better fit to available observations. 
The best-fit model fits the data (BAO+H(z)+Union) with x 2 = 317.64 (for 314 degrees 
of freedom) which means that the model can be ruled out at only the 56.8% level. The 
best-fit model has H 5h ~ 8.5 km s _1 Mpc -1 and a ~ 1.3 Gpc. We note that this 
large scale flow is not compatible with the original Hubble Bubble from Ref. [13] where 
H Sh ~ 4.5 km s _1 Mpc -1 and a ~ 0.1 Gpc. Indeed, as seen from Fig. [5] the Hubble 
Bubble can be ruled out at 2a confidence. However, when the Riess gold data set is 
used to constrain the model, the Hubble Bubble is only ruled out at the 1-sigma level. 
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Table 1. Expected values (means) for H , a, and S p based on the BAO+H(z)+Union 
data sets. 



model H a a 

[km s' 1 Mpc" 1 ] [Gpc] 



cosmic void 
cosmic flow 



67.62 
69.24. 



+0.49 
-0.48 
-0.66 
-0.62 



2.96 



+0.19 
-0.45 



4.05+i; 



44 

i 4fi +0.55 niQ+0.02 
i -^°-0.26 u - lo -0.01 



Table 2. Parameters and the x 2 of the best-fit models based on the BAO+H(z)+Union 
data sets. 



model 


[km s" 1 Mpc -1 ] 


a 
[Gpc] 


6 


X 2 


X 2 /dof 


cosmic void 
cosmic flow 


67.603 
69.209 


3.119 
1.327 


4.501 
0.123 


336.55 
317.64 


1.07 
1.01 



2000 



1500 



o 



^ 1000 



< 

Q 



500 



/ ^ 

/ < 


0.02 

-0.02 
-0.04 
-0.06 
-0.08 

-0.1 






..„„„„ 




bf 1 


bf2 
ACDM 

1 
















~--^^ 








D 


0.2 


0.4 


0.6 


0.8 



0.5 



1.5 



2.5 



Figure 6. The angular diameter distance as a function of redshift in the ACDM 
model and the best-fit inhomogeneous models (bf 1 - best-fit cosmic void model, bf 2 
- best-fit cosmic flow model). The inset presents the absolute distance difference, i.e. 
{D — -Dacdm)/-D (where D is the angular distance either in bf 1 or bf 2 model). 



6. Future observational constraints 



The previous section showed that very large regions of parameter space for 
inhomogeneous models can be ruled out by combinations of existing constraints. 
However, a large scale void model or an ultra-large cosmic flow remain as possible 
alternatives to dark energy with respect to explaining the apparent acceleration of 
the Universe. In our analysis we did not consider constraints from the CMB. This 
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Figure 7. The residual Hubble diagrams for the best-fit void model (bf 1), best-fit 
cosmic flow model (bf 2), and the ACDM model. 



is because observations at high redshift (z ~^> 1) are of limited significance when used 
to constrain the properties of the local-scale inhomogeneities. Owing to the flexibility 
of inhomogeneous models we can always assume that our void is surrounded by an 
overdense or underdense region - there is not need to assume that the Universe consists 
only of the void of Gpc-radius and everywhere else is homogeneous. Thus, by varying the 
density of such a ring or its size, every model studied in this paper can be made to fit the 
distance to the last scattering surface (which is given by the position of the first of the 
CMB peak). Thus, the high redshift data does not significantly constrain the properties 
of a local inhomogeneity (see Appendix for more details). To tighten the constraints on 
the properties of the local Gpc-scale inhomogeneity we need more precise data sets at low 
redshift which can be directly compared with constraints from existing data sets such as 
Type-la supernovae. Among these observations, the most valuable are those of distance, 
of the Hubble parameter, or measurements like baryonic acoustic oscillations and the 
Alcock-Paczyhski tests which depend on both the Hubble parameter and distance. For 
example, precise measurements of a maximum in the angular diameter distance distance 
would place very tight constraints [211 [22]. This is depicted in Fig. [6] which shows the 
angular distance as a function of redshift in the ACDM model and in the best fitted 
models for both families of models studied in this paper. The angular distance to 
z > 1 in inhomogeneous void models is distinguishable from the corresponding distance 
in homogeneous models, at a level larger then 10%. In addition and the position of 
maximum for model bf 1 is shifted by Az « 0.5 relative to the ACDM model. However, 
current estimates of the angular distance from Snla are predominantly at z < 1, and 
are not sufficiently precise to discriminate between models on this basis. Other state 
of the art the measurements of the distance based on the Sunyaev-Zel'dovich effect are 
also too imprecise [23]. Another possible test that has been proposed in the literature is 
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based on the time drift of cosmological redshift, i.e. a measurement of z [24]. However, 
this requires a very precise measurement since the typical amplitude of the effect is 
8z ~ 1CT 10 on a time scale of 10 years [24]. Another possible test was proposed in [2-5J , 
which relies on studying the consistency relation between the distance on the null cone 
and the Hubble parameter. Such relation does not hold if the geometry of the Universe 
differs from the FLRW geometry. Of more immediate utility are methods based on 
spectral distortions of the CMB power spectrum [26]. These have already been used 
to put an upper bound of ~ 1 Gpc for the radius of a local Gpc-scale void [27]. In 
addition, a study of the kinematic Sunyaev-Zel'dovich effect has recently been used 
by Garcia-Bellido and Haugb0lle [7] to estimate the properties of a local void. This 
analysis shows that the void permitted by the existing data cannot be larger than ~ 1.5 
Gpc. We recall from Sec. 15.11 that a void with a = 1 Gpc and 5 P = 0.95 fits the data 



70 



with x — 347.27, which means that that the model can be rejected at only the 90. 1 
level, which is insufficient to rule out this possibility. Moreover, cosmic-flow models, 
although not tested with these type of observations (kSZ effect and spectral distortion), 
require an inhomogeneity of radius also around 1 Gpc. Therefore, there is a need for 
future observational constraints to either confirm or rule out the existence of these large 
inhomogeneities. 

Of considerable promise for the future are very precise measurements of supernova 
in the redshift range of 0.1-0.4 [8]. In this redshift range (see Fig. [7]) measurements of 
over 2000 supernova (for example by the Joint Dark Energy Mission) will be sufficient 
to reassure existence of such a void. However, in the near future the most promising 
candidate for testing inhomogeneous cosmological models are BAO measurements. 
Current measurements of BAO from LRG SDSS (for an up to date status see [28] ) 
already add additional constrains (Figs. [3] and [5]). Moreover, the WiggleZ project |§| will 
soon measure the dilation distance with accuracy of 2.5% at z = 0.75. In the next 
subsection we show that this forthcoming observation will be sufficient to significantly 
constrain the inhomogeneous cosmological models. 

6.1. BAO constraints 

Before proceeding with further analysis, we asses the importance of analysing BAO data 
using inhomogeneous, rather than FLRW models. For example, when analysing the data 
in pjj] the redshifts of LRG galaxies were translated into comoving coordinates using 
the flat FLRW model with Q m = 0.3, Q A = 0.7, H = 70 km s" 1 Mpc" 1 . However, if the 
ACDM model is not a good background model of the Universe then such a procedure 
leads to inaccurate results. The error which arises from this conversion can be estimated 
from the inset in Fig. [H As seen at z = 0.35, the error in the distance conversion is 3% 
and 0.3% respectively in cases where the best-fit model from the first (bf 1) and second 
(bf 2) family of models studied in this paper were chosen as the background rather 
than the ACDM model. These results indicate that while the effect is relatively small, 

§ http://wigglez.swin.edu.au 
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we should be aware that the value Dy at z = 0.35 estimated in [19] when using the 
inhomogeneous model could be different than Dy = 1370 ± 64 Mpc. On the other hand 
we find that even if this value was change by 5% the results presented in the previous 
section would still hold. This is because the, as seen from Figs. [3] and [5] the observations 
of BAO at z = 0.35 do not put tight constraints on the inhomogeneous models relative to 
data. However at higher redshift (z~ 1) the redshift-distance conversion error increases 
to around 10% for the cosmic void model. Thus, when using forthcoming BAO data to 
constrain inhomogeneous models one must perform with the whole analyses within the 
inhomogeneous framework. 

In this section we will assume for the sake of argument that such an analysis has 
been performed. We will also assume that the Dy is equal to the value obtained within 
the ACDM model (Q m = 0.3, Q A = 0.7, H = 70 km s" 1 Mpc" 1 ). As in Sec. E 
we set a G [0,6]Gpc, 5 P G [0,6], and 5h G [0,0.4]. However, we tighten H range 
to H G [68,72], because future observational experiments like WiggleZ project will 
measure the Hubble parameter with 4% precision. 

The potential constraints from this BAO experiment around the ACDM model are 
presented in the left panels of Fig. El As measurement of the BAO scale at z = 0.75 
with a precision of 2.5% would rule out most of the parameter space considered for the 
cosmic depression family of models. A measurement at z = 0.95 with 2.5% accuracy 
would be sufficient to rule out all remaining possibilities in this case at a significance of 
3a. 

The second family of models considered in this paper are not tightly constrained 
by a single BAO measurement at z — 0.75. As seen in lower right panel of Fig. [8] even 
measurements up to redshift 1.75 are not able to rule out this type of model. To rule 
out all remaining possibilities at a significance of 3<r measurements should be extended 
beyond z = 2. 

7. Conclusions 

Homogeneous cosmological models which include a dark energy component and which 
are built on the Copernican Principle have been spectacularly successful in describing a 
suite of cosmological observations. However, alternatives which replace dark energy 
with an gigaparsec-scale inhomogeneity have also been postulated as explanations 
for the apparent acceleration of the Universe. In this paper we have shown that 
most inhomogeneous anti-Copernican models can already be ruled out by combining 
current cosmological observations (SNIa, H(z), and BAO). Our model were based on 
a simple parametrization of an inhomogeneity. We assumed that the inhomogeneity is 
modeled with a Gaussian-like profile, although other profiles are also possible. For 
example a similar analysis was recently proposed by Garcia-Bellido and Haugb0lle 
[30] were the inhomogeneity was modeled using a hyperbolic tangent function. Their 
results also support our conclusion that there is still a range of parameter space that 
remains consistent with current constraints. As a result the success of Concordance 
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Figure 8. Observational constraints on the parameters of the first family of models, 
shown as contour plots of x 2 f° r the model given the data, after marginalizing over 
Hq. The confidence levels are set assuming that the value of Dy is the same as in 
the ACDM model and that it is measured with 2.5% precision. Upper Left: expected 
constraints from the BAO at z=0.75 for the first family of models (cosmic void). Upper 
Right: expected constraints from the BAO at z=0.95 for the first family of models 
(cosmic void). Lower Left: expected constraints from the BAO at z=0.75 for the 
second family of models (cosmic flow) . Lower Right: expected constraints from the 8 
BAO measurements at z=0. 35,0. 55,0.75,. ..,1.75. 



models cannot be considered as a validation of the Copernican Principle. Forthcoming 
experiments will be capable of constraining all possible alternatives and will therefore 
either validate the Copernican Principle or determine the existence of a local Gpc-scale 
inhomogeneity. 
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Appendix A. Fitting the CMB data with inhomogeneous models 

In the standard approach, the CMB temperature fluctuations are analyzed by solving the 
Boltzmann equation within linear perturbation around the homogeneous and isotropic 
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FLRW model [HH E2JJJJ. Within an inhomogeneous background one can proceed with a 
similar analysis - employing the LT model instead of the FLRW model. Alternatively, 
if it is assumed that the early Universe (before and up to the last scattering instant) is 
well described by the FLRW model then the CMB power spectrum can be parametrized 

by m 



' r, 



l a (m - (j) m ) 



where l±, I2, h is a position of the first, second and third peak, and 
R 



In 



71- 



(A.l) 



(A.2) 



where R is a comoving distance to the last scattering surface and r s is a size of the sound 
horizon at the last scattering instant. The sound horizon depends on Q b , Q m , fl 7 , and h 
(where fl&, Q m , f2 7 are ratios of baryon, mater and radiation energy density respectively 
to the critical energy density evaluated at the current instant and h is equal to H/100 
km s _1 Mpc -1 ). The size of the sound horizon can be calculated accordingly to formulae 
given in [34]. The function <p m depends on f^, Q m , h, n s , and the energy density of 
dark energy (n s is the spectral index). The exact form of 4> m is given in [33]. The 
Vt m and h are given by a ratio p/p cr and H/(100 km s~ x Mpc -1 ) at the last scattering 
surface (evaluated for the current instant). 

In this appendix we will present a fit to the CMB data using the above formula. 
Let us consider four models 



model Al: 

p(t , r) = p b 1 + 1.45 - 1.45 exp 
t B = 0, tt b = 0.0445, n s = 1, 
where £ = r/Gpc. 

model Al+ring: 



Pb 



1 + 1.45- 1.45 exp 



0.08, 



n x 



p{to,r) 

t B = o, 

model A2 

p(t ,r) = p b 3.33 + 1.4 exp 
t B = 0, Q b = 0.0445, n s 
model A2+ring: 



0.963. 



0.75 2 
1. 



P(t ,r) 
t B = 0, 



Pb 

n b 



3.33 + 1.4 exp 



0.75 2 

0.07, n, = 0.963. 



c 2 



0.75- 



0.75- 



1.75 exp {-(y#y 



+ 0.08 exp 



\ V 0.965 I 



The density distributions of these models are presented in Fig. IA11 The fit to 



positions of the CMB peaks is presented in Table Appendix A As can be seen, model 



This 



approach 



is 



implemented 



such 



codes 



like 



(http://www.cfa.harvard.edu/~mzaldarr/CMBFAST/cmbfast.html), 
(http://www.camb.info/), or CMBEASY (http://www.cmbeasy.org/) 



CMBFAST 
CAMB 
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Table Al. CMB fit 

Model First peak Second peak Third peak 



Al 
A2 

Al + ring 
A2 + ring 



177.34 


414.08 


627.19 


216.59 


503.37 


754.28 


220.46 


530.38 


800.19 


220.43 


530.51 


779.76 



WMAP [351 H] 220.8 ±0.7 530.9 ± 3.8 700-1000 



Al does not fit the observed CMB power spectrum. However, we can modify the density 
distribution is such a way (by adding an underdense ring) that the CMB peaks can be 
reproduced - model Al+ring. In models where density increases up to some distance 
from the origin (local void type of models) a satisfactory fit to the first peak of CMB 
power spectrum (i.e. distance to the last scattering instant) can be obtained if the mass 
of the universe is decreased. This decrease can be obtained either by having a local void 
of a larger radius (as in [6] where R ~ 2.5 Gpc), or by having an addition underdense 
ring between local void and the last scattering surface. 

On the other hand, it is not only the mass that is important. To show this let us 
consider model A2. In model A2 the density decreases up to some distance from the 
origin. In this configuration we can also obtain a good fit to the position of the first 
peak if we increase the mass of the universe. 

The positions of other peaks strongly depend on Qf,. To obtain a good fit to the 
positions of other peaks in our model we need to increase the value of fib beyond the value 
that is consistent with observation of light elements in the local Universe (fib = 0.0445 
|36j) to Qb = 0.08. However, due to the lower expansion rate at large distances, the 
physical baryon density flbh 2 remains is close to the observed value, i.e. 0.0198 and 
0.0212 for models Al+ring and A2+ring respectively. 

The above results suggest that almost every model can be modified in such a way 
that a good fit to CMB power spectrum data can be obtained. Models Al and A2 are 
very different and yet after some modifications they are both able to fit the CMB data. 
Thus, the CMB data does not strongly constrain the properties of the local Gpc void. 
Each model for the local scale void considered in this paper can be modify in such a 
way (by adding one or two rings between a local void and the last scattering surface) 
so that the CMB power spectrum is recovered. 
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